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Abstract 

Effects of CP violation on the supersymmetric electro-weak correction to the 
anomalous magnetic moment of the muon are investigated with the most general 
allowed set of CP violating phases in MSSM. The analysis includes contributions 
from the chargino and the neutralino exchanges to the muon anomaly. The su- 
persymmetric contributions depend only on specific combinations of CP phases. 
The independent set of such phases is classified. We analyse the effects of the 
phases under the EDM constraints and show that large CP violating phases can 
drastically affect the magnitude of the supersymmetric electro-weak contribution 
to and may even affect its overall sign. 



1 Introduction 



As is well known supersymmetric theories contain many new sources of CP viola- 
tion which mostly arise from the phases of the soft SUSY breaking parameters and 
that such phases contribute to the electric dipole moments (EMDs) of the electron 
and of the neutron. Experimentally the electron and the neutron EDMs have very 
strict limits, i.e., for the neutron the limit is[|ll 



\dn\ < 6.3 X IQ-^^ecm (1) 
and for the electron the limit is|0] 

141 < 4.3 X 10-2^ecm (2) 

and these limits impose stringent constrains on particle physics models. In SUSY/ string 
models one normally expects CP violating phases 0(1) and phases of this size typ- 
ically lead to EDM predictions in such models already in excess of the current 
experimental limits. Of the possible remedies to this problem the conventional 
approach has been to assume that the phases are small |0, ^, typically, 0(10^^^), 
which, however, constitutes a fine tuning. Another possibility suggested is to as- 
sume that the SUSY spectrum is heavy in the several TeV region|p. Generally, a 
heavy spectrum may constitutes fine tuning]^ except in certain limited domains 
of the parameter space |[. Further, such a heavy spectrum may he outside the 
reach of even the Large Hadron Collider (LHC) and thus a disappointing scenario 
from the point of view of particle physicists. A third more encouraging possibility 
is that the large phases could indeed be there, but one escapes the experimen- 
tal EDM constraints because of cancellations among the various contributions to 
the EDMs. This possibility was proposed in Ref.[^ and there have been further 



verification and developments ||10|, |12|, |13|, |14| and applications such as in dark 
matter [ll5|, llB, |18|, [T^, E0[' energy processes pi], E2|, |2^, E4i, and on other 



SUSY phenomena H, |6], |3, |2|, |2|, |30 



The cancellation mechanism opens a new window on the SUSY parameter 
space where large CP phases along with a light SUSY spectrum can co-exist. 
Thus significant effects on SUSY phenomena can result. One of the quantities 
affected by CP phases is = {g^j, — 2)/2, where — 2 is the anomalous magnetic 
moment of the muon. This quantity is of considerable current interest since the new 
Brookhaven experiment will measure to an accuracy of better than a factor 
of 20 p^. Further, recently there has been considerable progress in reducing the 
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hadronic error |36|. With the reduced hadronic error the new experiment 

will test the Standard Model (SM) electro-weak correction]^ which including the 
two loop SM corrections stands atp8[ 



a^^ = 15.1 X 10-1° (3) 



It turn out that the supersymmetric electro-weak corrections to can be quite 
large and these supersymmetric effects on have been investigated for many 
years|3^, ^ However, the CP violating supersymmetric electro-weak 

effects on have been ignored for the reason that small CP phases or large CP 
phases with a heavy spectrum lead only to negligible effects on a^. 

With the cancellation mechanism the possibility of large CP phases along with 
a light spectrum arises and such a situation can lead to very significant effects 
on a^. Indeed in a recent work|^ the effects of CP phases on in the context 
of the minimal supergravity model (mSUGRA) were analysed and it was shown 
that CP violating phases can produce significant effects on a^. In the absence 
of CP violating phases the soft SUSY breaking parameters at the GUT scale 
in mSUGRA|Q consist of the universal scalar mass mg, the universal gaugino 
mass mi, the universal trilinear coupling Aq and tan/5 =< H2 > / < Hi > 
where H2 gives mass to the up quark and Hi gives mass to the down quark. 
More generally the soft SUSY breaking parameters as well as the Higgs VEVs 
are complex and have phases. However, by a redefinition of fields it is easily seen 
that there are only two CP violating phases in mSUGRA. These can be chosen 
to be the phase of Aq and the phase of /xq where /iq appears in the Higgs mixing 
term, i.e., in the term ^qHiH2 in the superpotential. In this paper we extend our 
analysis of the effects of CP violating phases on to supergravity models with 



non- universalities ||45|, |47| , ^ pO| , |5T| , |52| and to the Minimal Supersymmetric 
Standard Model (MSSM) which has many more CP violating phases. The existence 
of a larger set of CP phases widens the region of the parameter space where 
cancellations can occur. The purpose of this paper is to derive the general one 
loop supersymmetric correction to with the most general set of CP violating 
phases allowed in MSSM and determine the numerical effects of these CP violating 
phases on under the experimental constraints on the electron and on the neutron 
EDM. 

The outline of the rest of the paper is as follows: In Sec. 2 we derive the general 
one loop formula for aj for the case of a fermion / interacting with a fermion and 
a scalar in the presence of CP violating phases and without any approximation 
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on the relative size of the external and internal particle masses in the loop. In 
Sec. 3 we apply this formula for the computation of the chargino and neutralino 
exchange contributions to for the most general allowed set of CP violating 
phases in this sector. In Sec. 4 and in Appendix A wc study the combination of 
CP phases that enter and compare them with the corresponding combinations 
that arise in the expressions for the electron and the neutron EDMs. In Sec. 5 and 
in Appendix B the supersymmetric limit of our result is given and it is explicitly 
shown how the one loop Standard Model contribution to including the one loop 
QED correction, i.e. aem/27r, is cancelled by the supersymmetric contribution. In 
Sec. 6 we give a discussion of the satisfaction of the EDM constraints. In Sec. 7 
we give an analysis of CP violating effects on a^. Conclusions are given in Sec. 8. 
Appendix C is devoted to a discussion of the vanishing CP violating phases and a 
comparison of our results with previous analyses. 

2 CP Effects on - 2 in MSSM 

We give here the general analysis for the CP effects on — 2 of a fermion. In 
general for the interaction of a fermion ipj oi mass m/ interacting with a fermion 
■01 of mass rrii and a scalar 0^ of mass mfe, the vertex interaction has the general 
form 

- Ant = YlilJfiKik — + — Y^)'4^i(t>k + H.c. (4) 

ik 

This interaction violates CP invariance iff Im{KikL*,^) ^ 0. The one loop contri- 
bution to a/ is given by 

o/ = Oy- + (5) 

where aj. and arise from Fig. 1(a) and Fig. 1(b) respectively. Oy- is a sum of 
two terms: = a^^ + of where 

=ll^^<^^^^k)h{%% (6) 



oTT'^mi mi mi 



and 



and where 



rl rl-x ^ 

hia, 13)^—1 dx dz — — (7) 

^ ^ Jo Jo az^ + {1 - a - /3)z + /3 



- E T^i\K^k\' + \L^km^, ^) (8) 

■' ^ loTr'^mf mf mf 
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and 



r-l r-l-x ^2 — ^ 

hia, p) = dx dz ^ ^^ , , ^ (9) 



q;2;2 + (1 - a - + 
Similarly, consists of two terms: aj — of- + of where 



2 ™2 



and 



and where 



and 



(10) 

^-'°-^)=/o'''^r'''- „.^ + (/>-"a-l). + l 
«? = - S + |i..R^4(5, 3) (12) 



/■I r\-x ^2 _ ^ 

74(0;, dx dz — — r — (13) 

^ ' Jo Jo az^ + {f3-a-l)z + l ^ ' 

In the above we have given the exact expressions for the integrals Ji —I4 rather than 
their approximate forms in the limit when one neglects terms of size m'j relative 
to ml and m| which allows one to write simple closed form expressions for them. 
We will see that the general expressions are needed to discuss the supersymmetric 
limit of our results which provides the absolute check on our normalizations. 

3 On — 2 with CP Violating Phases 

We apply now the above relations for the computation of the chargino and the neu- 
tralino exchange contributions. We consider the chargino exchange contributions 
first. The CP violating phases enter here via the chargino mass matrix defined by 

where xi ^-^id X2 are phases of the Higgs VEVs, i.e., < Hi >— \ < Hi > \e^^^ 
(i=l,2). The matrix of Eq.(14) can be diagonahzed by the biunitary transformation 
U*McV~^ — diag{m.^+,rh^+) where U and V are unitary matrices. By looking at 
the muon-chargino-sneutrino interaction one can identify Ki and Li and one finds 

< = < + < (15) 
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where aj^ and a^^ are given below. We exhibit these only in the limit where /a (a, (3) 
and 74(0;, /3) have their first arguments zero and one may write 

h{0,x) = ~Fs{x), 140, x) = -^F,{x) (16) 



where 



F^ix) = -. TT^fSx^ - 4x + 1 - 2x'^lnx) (17) 

(x — 1)'* 

F^x) = - — ^—-^{2x^ + - 6x + 1 - Gx'^lnx). (18) 

(a; — 1) 



In the above approximation we have 



and 



ET?^(l'^M^i^r + l^ar)F4(x7f)- (20) 



where 



I — 1 vT V. 
^ ^ "^M ^-ixi (21) 

" V2Mw COS P ^ ' 

Next we discuss the ncutralino exchange contribution to a^. CP violating 
effects here arc all contained in the neutralino and smuon mass matrices. For the 
neutralino mass matrix the CP violating phases enter as below 

( |mi|e^^i -M^ sin cos /9e-*» sin ^vi^ sin /?e-^>^2 \ 

|m2|e*«2 M^cos^H^cos/?e-^» -M^cos^vk sin/^e-^'^^ 

sin cos /5e-^» cos cos /^e'^^^ q -|/x|e*^'^ 

V sin sin /?e-^^i -M^ cos sin /Je'^^^ -l/xle^^" / 

(22) 

The neutralino mass matrix M^o is a complex non hermitian and symmetric matrix 

and can be diagonalized using a unitary matrix X such that X-^M^oX=diag(mj^o, m^^o, m^o, m^o). 

Since the loop correction involving the neutralino exchange also involves the smuon 

exchange (see Fig. la) the CP phases in the smuon (massy also enter the analysis. 

The smuon [massY matrix is given by 

^2^( Ml,, m^(A;mo-/xtan/3e*(>^i+>^^))\ 

This matrix is hermitian and can be diagonalized by the unitary transformation 

L'tM|L> = diag(M|i,iW|2) (24) 
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The neutralino exchange contribution to is given by 

«^ = af + af (25) 



where 



and 



a 



11 TJii^aEM A 1 ^ , I. ,ml Mj- 



.-,,35,., i:g^«e(4-)«4.4) (26) 



rilj = -(-^[tan^H^Xi, + Xs,]/^^, - k^X:,,DI^) 

(V^tan 9wXijD2k + K^X^jDik) (27) 

and a^^ is given by 

,12 _ rnlaEM 'J^ 1 r ( ^ooN 



where 



2 

^ I Y |2 

'^^^ 2M2,cos2/?l 

+^tan2%|Xi,f(|L'i,|2 + 4|L>2fer) + ^|X2,f|i^ifer 

+ tan^w^|Difcpi?e(Xi,X2*,) 

-Re{e-'^^XsjX*.D,kD;,) (29) 



M^^ cos P 

If one ignores the muon mass with respect to the other masses involved in the 
problem, the form factors Ii{a,P) and l2{a,P) become 

h{0,x) = ^Fi(x),/2(0,x) = ^F^ix) (30) 



where 



and 



Fi{x)= ^ {1 - x"^ + 2xlnx) (31) 
[x 1) 

F2{x) = y — + Gx^ - 3x - 2 - GxZnx). (32) 
[x — 1) 
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4 The number of independent linear combina- 
tions of phases that enter 

Not all the phases that enter in the chargino, neutralino and smuon mass matrices 
are independent. We discuss here the set of independent phases that enter a^. 
We consider the chargino contribution to first. Here the matrix elements of U 
and V as defined in the paragraph following Eq.(14) along with k,^ as defined by 
Eq.(21) carry the phases ^2, 0^, xi and X2- By introducing the transformation 
Mc = BrM'(jB{ and choosing Br = diag{e^^^ , e''^^) and Bl = o?ia^(l, e*(»+«2)) 
we can rotate the phases so that is given by 

M> -( l"^2| ^/2mw fivn [3 \ , . 

The matrix Mq can be diagonahzed by the biunitary transformation C/^M^C/^^diag 
(m^+, m^+). It is clear that the matrix elements of C/^ and Ur are functions only of 
the combination 6 = Ofj, + ^2 + Xi + X2- We also have U*McV~^ — diag{m^+,m^+) 
where U — (BrUr)^, and y^—{B]JJL)^. By inserting the new forms of U and V 
in the chargino contribution one finds (as shown in Appendix A) that a?j} and a^^ 
depend on only one combination, i.e., ^ = + ^2 + Xi + X2- 

Now we turn to the neutralino contribution, the phases that enter here are 6^, 
(^A^-i ^2, ^1, Xi and X2 and they are carried by the matrix elements of X, D^j^ and 
the phase of k^. Next we make the transformation M^o—P^o M'^o Px° where 

P^o = diag{e''-^,e''-^,e-'^'-^+^'\ e-*(^+>^2)) (34) 
After the transformation the matrix M^o takes the form 



/ \rhi\ —Mz sin 6w cos (3 sin 6w sin j3e~'^~ \ 



\fri2\ Mz cos 9 w cos (5 —M^ cos 6*14/ sin /5 

— M2 sin6'vi/ cos/5 cos^viz cos/3e*^ — |/x|e*^' 

V sin 6'vi/ sin /3e~*^ — cos sin /3 — |/x|e'^' 

(35) 

where 9' = + + xi + X2, and = (^1 - ^2)- The matrix M^o can be 
diagonalized by the transformation y-^M^o5^=diag(m^o, m^o, m^o, m^o) where Y 
is a function only of 9' and Thus the complex non hermitian and symmetric 

matrix M^o can be diagonalized using a unitary matrix X — P^oY such that 
X'^M^oX — diag(m^o, m^o, m^o, m^o). As shown in Appendix A by applying 
the above transformations to each term of t;^^ and X^^ one finds that the phase 
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combinations that enter here are Q\ and aA^, + + Xi + X2 from which we 

can construct the three combinations G + + Xi + Xi-, 'C2 + 6'^^ + Xi + Xi and 

By defining Q\ = ^^/^ + Xi + X2 one finds that the dependence on phases 
from both the chargino and the neutrahno exchanges consists only of the three 
combinations aA^, +6*1, G + ^1 ^^'^ '62 + ^1- One may compare these combinations 
with those that appear in the supersymmetric contribution to the electron and 
the neutron EDMs. In the analysis of Ref . |TD[ we found that the electron and the 



neutron EDMs depend on the following combinations: + ^^i, [i = 1,2,3) and 
c^Ak + ^1 with k = u,d,t,b,c, s; I. We note that even though and the EDMs 
are very different physical quantities the linear combination of phases that enter 
in them are similar. In fact the phases that enter subset of phases that 

enter in the supersymmetric contributions to the EDMs. 

5 The Supersymmetric Limit 

To check the absolute normalization of our results we discuss now their supersym- 



metric limit. In Ref. P3| the supersymetric contributions to from the chargino 



sector in the supersymmetric limit were computed by going to the limit such that 

U*McV-^ = diag{Mw, Mw) (36) 

In this limit it was shown that the contributions from this sector was precisely 
negative of the contribution from the W exchange |P7||. The analysis of Ref.p3[ 



was carried out in the framework of mSUGRA with two CP violating phases. For 
the MSSM case being discussed here with many CP phases the structure of U and 
V matrices in susy limit will be modified so that 

f-^(_\ :-:::). ''-^^(i :;x) pt) 

However, taking the phases Xi and X2 into account we find exactly the same result 
as in Ref.[|^ due to the appearance of the phase. Thus the sum of the W 



exchange contribution and of the chargino exchange contributions cancel in the 



supersymmetric limit. Similarly it was shown in Ref. by making a unitary 
transformation that the neutralino mass matrix in the supersymmetric limit can 
be written in the form 
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M^oX = diag{0, 0, Mz, Mz) (38) 

where the eigen-values are positive definite. It was then shown that the last two 
eigen-modes give a contribution which is negative of the contribution from the 
Z exchange in the Standard Model[0. In the case of MSSM we are discussing 
here the structure of the diagonalizing matrix X is now changed because of the 
Xi and X2 phases (see Appendix B). However, the final result we arrived at in 



Ref .11431 still holds due to the appearance of the phase once again. Thus the 
sum of the Z exchange and of the heavy neutralino exchanges exactly cancel in the 
supersymmetric limit. 

We now turn to the supersymmetric limit of the contribution from the first two 
massless eigen states of the neutralino mass matrix. A direct sum over the first 
two eigen-modes for the case of Eq.(38) in the supersymmetric limit gives 

ar^(^ero - modes) = (39) 

Thus the sum of the Standard Model contributions to from the photon at one 
loop and form the Z and the W exchanges at one loop 1 37 1 is cancelled by the 
supersymmetric contributions from the neutralino and the chargino exchanges in 
MSSM at one loop in the supersymmetric limit, i.e., in the supersymmetric limit 
one has 



^MSSM ^ Q ^4Q) 



This result is consistent with the expectation on general grounds The 
details of the derivation of Eq.(39) are given in Appendix B. 



6 Satisfaction of EDM Constraints 

Before proceeding to discuss the CP effects on we describe briefly the EDM 
constraints on the CP violating phases. As is well known for the case of the neutron 
EDM there are three operators that contribute to the neutron EDM, namely, the 
electric dipole moment operator, the color dipole moment operator and the purely 
gluonic dimension six operator. Both the electric and color operators have three 
components each from the chargino, neutralino and gluino contributions. For the 
electron case we have only the electric dipole moment operator which has only two 
components, the chargino and the neutralino ones. Recently, it has been pointed 
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out that in addition to the above contributions certain two loop graphs may also 
contribute significantly in some regions of the parameter space ||5^. In our analysis 
here we include the effects of these contributions as well. However, the effect of 
these terms is found to be generally very small compared to the other contributions 
in most of the parameter space we consider. Satisfaction of the EDM constraints 
can be achieved in a straightforward fashion using the cancellation mechanism. 



7 Analysis of CP Violating Effects 

In the above we have given the most general analysis of within the framework 
of MSSM with inclusion of CP violating phases. Our results limit to those of 
Refs. [p9t ^ in the limit when CP violating effects vanish (see Appendix C for 
details). For the case of the general analysis with phases in MSSM the number 
of parameters that enter along with the number of parameters that enter the 
EDM constraints which must be imposed on the CP violating phases is large. For 
the purpose of a numerical study of the CP violating effects on we shall confine 
ourselves to a more constrained set. Here we shall generate the masses of the 
sparticles at low energy starting with parameters at the GUT scale and evolve 
these downwards using renormalization group equations. At the GUT scale we 
shall use the parameters mo, mi/2, and Aq, and /i will be determined via radiative 
breaking of the electro- weak symmetry. We set Xi + X2 = and choose the phases 
that we vary to consist of 9^, Q-Aqi and {i = 1,2,3). The choice of the above 
constrained set is simply for the purpose of reducing the number of parameters for 
the numerical study. 

We begin our discussion of the numerical results by exhibiting the dependence 
of on the CP violating phases but without the imposition of the EDM con- 



straints. The dependence of on 9^ and was already studied in Ref.pS 
and we confine ourselves here to the dependence of on ,^1 and .^2- In Fig-2 we 
exhibit the dependence of on ^1 and in Fig. 3 the dependence of on ^2- From 
Figs. 2 and 3 we find that is significantly affected by the dependence of both ^1 
and .^2- However, a comparison of Fig. 2 and Fig.3 shows that the dependence of 

on ^2 is much stronger than on ^1. The reason behind this difference is easily 
understood. The relatively weaker dependence on the ^1 phase arises because this 
phase appears only in the neutralino contribution while the ^2 phase appears both 
in the neutralino and in the chargino contributions to a^. 

We discuss now the effects of CP violating phases on under the EDM con- 
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straints. In Table 1 we show four points which he on the curves of Fig. 2 and Fig. 3. 
As one can see the SUSY mass parameters mo and mi/2 are relatively small (i.e., 
mo,wii/2 << 1 TeV), the CP phases are large and there is compatability with 
experimental constraints on the electron EDM and on the neutron EDM as a con- 
sequence of the cancellation mechanism. One also finds on comparing with and 
without phases that the effects of CP violating phases on are very significant. 



Table 1: 









c/„( 10-26 ecm) 


4(10-2Vm) 


[afj,{phases)]{10 ^) 


k(o)](io-^) 


(1) 


2.35 


.4 


-3.08 


-0.86 


-4.8 


7.45 


(2) 


1.98 


0.4 


-0.34 


-1.67 


-7.8 


11.7 


(3) 


1.2 


-1.5 


1.87 


2.24 


-3.25 


5.6 


(4) 


2.7 


-0.4 


1.87 


-0.03 


-15.5 


3.15 



Table caption: Parameters other than those exhibited corresponding to the cases 
(l)-(4) are: (1) mo=70, mi/2=99, tan(3=3 , |v4o|=5.6, ^i=-l, 6=1-5, 6=0-62; 

(2) mo=80, mi/2=99, tanf3=5 , |Ao|=5.5, 6=-0-8, 6=1-5, 6=0-95; (3) mo=75, 
mi/2=132, tan/3=A , |Ao|=6.6, 6=-l, 6=1-78, 6=2-74; (4) mo=70, mi/2=99, 

tanP=6 , |74o|=3.2, 6 =0-63, 6=0-41, 6=0-47, where aU masses are in GeV units 

and all phases are in rad- 



In Fig. 4 we exhibit Ci^ clS 8b function of mi/2 where all points on these trajectories 
satisfy the experimental constraints on the electron EDM and on the neutron EDM 
by cancellation. One finds that the magnitude of the supersymmetric electro-weak 
contributions are comparable to and even larger than the Standard Model electro- 
weak contribution as given by Eq.(3)|]38|. 



8 Conclusions 

We have given in this paper a complete one loop analysis of the effects of CP 
violating phases on with the most general set of allowed phases in MSSM in 
this sector. We have checked the absolute normalization of our results exhibit- 
ing the complete cancellation of the supersymmetric result in the supersymmetric 
limit with the Standard Model result including the qed one loop correction to a^, 
i.e., aem/27r. A detailed numerical analysis of the CP violating effects on for 
the regions which satisfy the EDM constraints is also given. Computations of 
under the EDM constraints shows that the supersymmetric electro-weak effects 
can generate significant contributions to even with moderate values of tan/3. 
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i.e., tan/9 ~ 3 — 6, which can be comparable to the Standard Model electro-weak 
correction. Thus supersymmetric CP effects on are within the realm of observ- 
ability in the new Brookhaven §^^ — 2 experiment. 
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Appendix A: 

In this appendix wc give the explicit derivation of the linear combinations of phases 
on which depends. For the chargino contributions the phases arc contained 
in the quantities Re{K^U*2V*^) and \k^JJ*^\^ + \Vii\^. By using U = (BrUr)'^ 
and V = {BlUl)^ as defined in Sec. 4, where Ul and Ur are functions of the 
combination = 6 ^ + ^2 + Xi + X21 one finds 



and 



which leads to 



and 



Uh-e^'^'U^R^i (41) 

= Ulu (42) 

^,U*2V:, = \K,\U*R2iULu (43) 

i^>.U;2? + \Vn? = K?\U*R2i\'' + \ULii\^ (44) 



Eqs.(43) and (44) show that the chargino contribution to depends only on one 
combination of phases, i.e., ^ = + + Xi + X2- 

For the case of the neutralino contribution to a^, the phases are contained in 
the quantities Re{r]^j) and X^j. We first consider the quantity Re{rj^j). It consists 
of six terms in the product 

{cXijD2k + K^XsjDik) (45) 

where a, b and c are real numbers and independent of phases. The first term in 
the expansion of Eq.(45) is 

acX^ Dl^D2k = ±acXlj cos / sin ^/e'^^ (46) 

where +(— ) sign is for /c = 1(2) and where the following definitions are used 

. 2m^[|"^o^Mp + lyU^Mp - 2|moA^/ii?/,| cosa]^/2 ( 
tan 20 f = (47) 



12 



Here = t&n.pe'^^^^'^^^^ and a = aA^^ + 0/^ + Xi + X2- The phase /?/ is defined 
such that ^ 

^"^^/= [A2 + ^2]l/2 (48) 

and 

'^f^f^WTW^^ (49) 

where A is defined by 

A = Imo^^l cos - l/^-R* | cos(^^ + Xi + X2) (50) 
and B is defined by 

B = \moA^\ shiaA^ + sin(6'^ + Xi + X2) (51) 

By using Xij — Yxje"^^^!"^ where Y^j are functions only of Q' and one can 

write the first term as given by Eq.(46) as follows 

acXlpX^D^k = acYl^h{a)e-'^^^-^^^ (52) 

where fkip) are real functions of a. By using the definition of /?/ as given by 
Eqs.(48) and (49) and by taking the real part of Eq.(52) we find that the right 
hand side of Eq.(52) contains the three combinations Q' . /S.(^/2 and a which come 
from the first part of the the right hand side of Eq.(52) and in addition it contains 
the following two combinations: ua^ — and + Xi + X2 + ^1 which come from 
the exponent. But the latter two combinations are linear combinations of the first 
three combinations. Thus the left hand side of Eq.(46) or Eq.(52) will depend only 
on the combinations 9\ and a. The same analysis can be applied to the other 

five terms and each one of them will give us the same three linear combinations. 

Next we consider X^^ . It consists of six terms and the quantities in them which 
contain phases are iXgjp, \Xij\'^{\Dik\^ + A\D2k\^), \X2j\^\Dik\^, \Dik\'^Re{XijX^j), 
Re{e-'^'X3jXljDikD*f.) and Re{e-'^'XsjX*jDikD*k). The first one of them, i.e. 
jXsjl^, can be written in terms of the Y matrix as [Isjp which depends only 
on the two combinations 9' and A^/2. The second expression can be written as 
\Yij\'^gk{a) where gk{oi) are real functions of a as defined after Eq.(47). So this 
term will depend on the three combinations 9\ A^/2 and a. The third expression 
is similar to the second one and will give the same combinations. The fourth 
expression can be written as hk{a)Re{YijY2j) where hk{a) are real functions of a. 
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So this term also depends on the same three combinations. The fifth expression 
can be written as 



Reie-'^'XsjXljD,kD*,) = ReiYsjY*jSkia)e'^^^-^f^) (53) 

where .s^.(«) are real functions of a. By treating the exponential term as we did in 
the first term of rj^j, it will give us two extra combinations besides the usual three. 
These are — Ci and 0/^ + xi + X2 + Ci which, however, are linear combinations 
of the usual three. Thus we end up here with the same three combinations. The 
sixth expression can be written as 

Re{e-'^^XsjX;^D,,D;,) = Re{YsjY;.s,{a)e'^'-^-f'f^) (54) 

from which we can identify the usual three combinations and in addition one has 
the following combinations in the exponent :aA^ — ^^-^ and 0^ + Xi + X2 + ^^-f^- 
These again are linear combination of the usual three and thus we end up with 
only three phases in the neutralino contribution, i.e., 9\ and a. 

Appendix B 

In this appendix we discuss the supersymmetric limit in the massless sector.. For 
this purpose we begin by exhibiting the unitary matrix X that diagonalizes the 
neutralino mass matrix in the supersymmetric limit such that the eigen- values are 
arranged so that 

X'^M^oX = diag(0, 0, Mz, Mz) (55) 
With the above ordering the unitary matrix X takes on the form 



/ 


a 




sin t)w 


■ sin 0w 






a tan 6w 


P tan 9\Y 


cos 9w 
V2 


■ cos e\v 










2^ 


2^ 




V 






2^ 


2^ 


1 



(56) 



where 



a = I — , 13 = —. — = (57) 



Using these results it is easily seen that the sum over the first two neutralino mass 
eigen- values gives 

a^i{zero - modes) = .T, ^EE ^^«.)(^) - (58) 



14 



where we set M^^ = and the factor H is defined by 



1 fl-X 



z 



H= dx dz- — (59) 

Jo Jo U — 1 



Thus in the supersymmetric hmit the entire supersymmetric contribution to a]^ 
from the masseles 
we need the sum 



from the masseless neutrahno states comes from a^^. To compute this contribution 



EEtttKMt^^T^) (60) 



where the sum over j runs only over the first two modes. In the supersymmetric 



o2 



hmit we set M^-^ = m^, M^o — > (j = l,2), and x^j = ^ ^ oo (j=l,2) and 

TTi^ 1 

-^hiXf^j^x^j) ^ -- (61) 



Now substitution of the exphcit form of the X matrix gives 



E E = 4^^^'^^ (62) 
j=i k=i 



Use of Eqs. (61) and (62) in Eq.(60) gives 



a 



a {zero — modes) = (63) 

^ 2n 

Thus we find that in the supersymmetric hmit the exchange of two massless neu- 
trahnos gives a one loop contribution to which is exactly negative of the photonic 
one loop contribution. Thus in the supersymmetric limit the sum of the one loop 
contributions of the zero modes of the theory cancel. The cancellation provides an 
absolute check on the normalization of our supersymmetric result in this sector. 
Appendix C 

In this section we consider the limit of vanishing CP violating phases and compare 
our results with those of previous works. We first compare our results with those 



of Ref.||39|. We consider the chargino contribution first. Using Eq.(2.8) of Ref.[]39 
and noting that the free part of the Lagrangian density for the complex scalar 
fields in that work is given by ^{d^z*d^z — m'^z*z), we find that our Ki and Lj are 



related to the and of Ref . ||3^ as follows: 



Ki,2u ~iV2At'- , Li,2. ^ -iV2A%^- (64) 
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Further, our form factors F^lx) and F4{x) are related to the form factors Fi and 
F2 of Ref.ra as follows: 



F.(x) 



-F.(x), FAx) 



FAx) 



(65) 



Defining 

gf = 2af, gf = 2af (66) 

we find that our Eq.(12) in the limit of vanishing CP violating phases is given in 
the notation of Ref. [B^] by 



9i 



w 



m. 



247r2 



4(1)2 .(a)2 



a=l,2 



rhl 



FAXa] 



(67) 



and similarly our Eq.(lO) in the same limit in the notation of Ref.[^ is given by 



m 



4772 ^ 



where 



a (a) 

^^F,{x^) 



a=l,2 



% a =1.2 

mi 



{61 



(69) 



Eqs. (67) and (68) agree precisely with Eqs.(2.6a) and (2.6b) of Ref. [p^ to leading 
order in /x^/m^ taking account of the typo in Eq.(2.6a) where A^^^ should read A^^^'^ 
and noting that J^]^~ is proportional to m^/M^ and thus does not contribute to 
leading order. 

We consider next the neutralino contribution. From the interaction Lagrangian 



Eq.(2.4) of Ref. p9[| we find the transition from our notation to that of Ref. |P9||as 
follows: 



where we identify O' to be 

O' - 



cos 6 sin 6 
— sin 6 cos 6 



(70) 



(71) 



Noting that our form factors Fi{x) and F2{x) are related to the form factors G2{x) 
and Gi{x) of Ref.pU| by 



Fi(x) = -G2(x), F2(x) = -2Gi(x) 



(72) 



and defining gf = 2a and g2 = 2a^ we find that our Eq.(8) gives precisely 



Eq.(2.10) of Ref.p9| taking account of the typos in Eq.(2.10a) in that 1/nk should 
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read l/j^l and G2{x2k) in the same equation should read Gi{x2k)- Further, our 



Eq.(6) agrees precisely with Eq.(2.12) of Ref.[^ 



Next we compare our results with those of Ref . |[40|| . For this purpose in the 
chargino sector we identify Wi and W2 states with the states W~ and H~ of 
Ref. PI in order to use Table 1 of Ref.|[4U||. With this identification in the limit 
of vanishing CP violating phases we find that in the chargino sector our matrices 
V and U are real and orthogonal and are related to the matrices Oi and O2 of 
Ref. HDI as follows 



^km ~^ Oimk, f^fcm ~^ ^2mfc (73) 

The analysis of Ref. computes only the contribution af^ of in their Eq.(5). 



Relating our -^3(77) to their Fsuirj) by F^ijf) 
can be written in the form 



-2Fsy{x), we find that our Eq.(19) 



2a 
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47r2sin2^^ ^ 



-p22kOlk^Fsu{yi. 



vk) 



\/2Mkmw cos (5 

which is exactly Eq.(5) of Ref.[^ on relating their smOu to our cos/9 by smOu 



(74) 



cos/?. The consistency of the analysis of Ref.|^0[ with our analysis, however, 
requires that the sign of the terms with My/ in the chargino mass matrix given by 



Eq. (3b) of Ref. |^ be reversed. 



To compare our results to those of Ref. |^ in the neutralino sector we note 
that our and H2 states are related to their and H^' states by = 
and H2 = H^' . In the limit of vanishing CP violating phases, our neutralino and 
smuon mass matrices become real, and the corresponding diagonalizing matrices X 
and D become orthogonal and can be identified with the real orthogonal matrices 
O and S of Ref. 0: 

X^O, D^S (75) 



The consistency of the analysis of Ref. |^ with our analysis, however, requires that 
the sign of the terms with Mz in the neutranlino mass matrix given by Eq.(3a) of 
Ref. 1^ be reversed. The analysis of Ref. calculated only the part aj} in their 
Eq.(6). To compare the result of aj} of our Eq.(25) with their Eq.(6) we first note 
that our Fi is related to their F by Fi{ri) = —Firf). Second we need to identify 
the fields Wi (i-1,2,3) in Eq.(6) of Ref.[|g in order to use Table 1 of Ref.[|| to 
write out in detail the interactions of Eq.(6). This identification is as follows: 



Wi = 5°, W2 = W3 = 



(76) 
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Further in Ref.|^ we identify L = 1 and R = 2 in their Eq.(6), and we need to 
complete their Table 1 since the term g^fiiH^s^) is missing in Table 1 and one 
needs it to expand out Eq.(6). Here we find that the entry for the magnitude for 
this coupling in their Table 1 should be the same as the magnitude for the coupling 
g{fijiH^sj^) listed in Table 1 (see Eqs. (5.1) and (5.4) of Ref.|]56|). Using the above 
correspondence we find that our result for 2aj} gotten from our Eq.(25) produces 
exactly Eq. (6) of Ref . [HDI in the limit of vanishing CP phases. 
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f x° f 

(a) 




(b) 



Figure 1: The one loop contribution to g"^ — 2 from (a) neutralino exchange, and 
(b) chargino exchange diagrams. 




^1 (rad.) 



Figure 2: Plot of a^*^ as a function of without the imposition of EDM con- 
straints. The values of the other parameters for the curves (l)-(4) correspond the 
cases (l)-(4) in Table 1. 



24 



50 — .- 

r.. 1. 



-50 







Figure 3: Plot of function of ^2 without the imposition of EDM con- 

straints. The values of the other parameters for the curves (l)-(4) correspond the 
cases (l)-(4) in Table 1. 
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Figure 4: Plot of a^^^^ as a function of mi where all points on the trajectories 
satisfy the current experimental constraints on the neutron and on the electron 
EDM. The curves labelled (l)-(3) are drawn for the following set of data: (1): 
|Ao|=6.5, ^^=2.92, a^,=-0.4, tan/? = 4, ^1 = 0, ^ = 0.2, ^3 = 0.065; (2): 
Uo|=5.4, 9^,= 3.006, «^o=-0.1, tan/? = 3.5, 6 = 0.105, 6 = 0.105, ^3 = 0.15; 
(3): |/lo|=2.9 , ^^=3.02 , a^o=0.5 , tan/? = 2.6, 6 = 0.19, ^2 = 0.19, {3 = 0.41. 
For all cases 50 < mo < 250 (GeV) and all phases are in rad. 
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